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Abstract 

For any integral operator K in the Schatten-von Neumann classes of 
compact operators and its approximated operator Kn (N > 1) obtained 
by using for example a quadrature or projection method, we show that 
the convergence of the approximate p-modified Fredholm determinants 
detN P (lN-\-zKN) to the p- modified Fredholm determinants &et v (I-u-\-zK) 
is uniform for all p > 1. As a result, we give the rate of convergences when 
evaluating at an eigenvalue or at an element of the resolvent set of K. 



1 Introduction 

Let K be an integral operator defined on a Hilbert space T~L and given by 

Ku(x) = / k(x,y)u(y)dy, x G X 
Jx 

where k(x,y) is such that for all x G X, 

k(x,.)u(.) G L X (X) and KueU. (1.1) 

Consider the following integral equation 

(I H + zK)u = f. (1.2) 

In this paper, we are not considering solving the above integral equation but 
rather the eigenvalue problem that is, when / = for all x G X. More precisely, 
we are interested in locating the eigenvalues of the operator K. For such a 
purpose, the approach used here is the numerical evaluation of the p-modified 
Fredholm determinants that is, det p (l^ + zK) for p > 1 [26l H5] . Indeed, 
the motivation behind this choice, comes from the fact that the the p-modified 
Fredholm determinants are entire functions whose reciprocal zeros are the eigen- 
values of the operator K with algebraic multiplicties equal to the order of the 
zeros [27] [I5j [26j [TTl 0. The first appearance of these functions, in particu- 
lar the 2-modified Fredholm determinants goes back to Hilbert [18] followed by 
Calerman [9], Smithies [27], Pelmelj [23] and others. Hilbert observed that in 
order to achieve a convergence of the determinant associated with an L 2 -kernel, 
it suffices to set k(x,x) to zero in the formula of the Fredholm determinant. On 
the other hand, the Fredholm determinant, as introduced by Fredholm [14], is 
an entire function of the spectral parameter which characterise the solvability 



of equation (1.2), under the assumption that the nonzero function / and the 
kernel k are continuous. Broadly speaking, the extension of Hilbert's theory 
to kernels which are not square-integrable has been developed by Gohberg and 
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Krein [16] and Dunford and Schwartz [13] . A particular attention for 4-modified 
Fredholm determinant was conducted by Brascamp [7]. Regarding the conver- 
gence issue, this also goes back to Hilbert. Indeed by applying the rectangular 
rule in equation (1.2), Hilbert showed that, for a continuous kernel function fc, 
the discrete version of the Fredholm determinant converges uniformly to the 
continuous one. Bornemann, in [8], extended Hilbert's result to any quadrature 
rule which converges for continuous functions. However in this paper, we prove 
the uniform convergence under weaker assumption of the kernel. Precisely, for 
kernel function k satisfying condition ( |1.1[ ) such that the corresponding integral 
operator K is in the Schatten-von Neumann classes of compact operators. Thus 
the desired convergence is obtained under the hypothesis that the set {K n } of 
numerical integral operators approximating K is collectively compact [1 . This 
is because if the approximation error (Ku—K n u) is uniform then the sequence of 
eigenvalues {A n } of K n converges uniformly to an eigenvalue A of K [T| |2] 13] [28] . 
The paper is organised as follows; in Section 2 we briefly review some proper- 
ties concerning the p- modified Fredholm determinants. In Section 3, we show 
that the p-modified finite dimensional determinants associated with the oper- 
ator K n converges uniformly to the p- modified Fredholm determinants. As a 
consequence, we give the rate of convergence when evaluating at an eigenvalue 
or at an element of the resolvent set of K. Finally in Section 4, we present some 
numerical results that demonstrate the analysis carried out in Section 3. 



2 The ^-modified Fredholm determinants 

We briefly recall some basics that we use in this paper. Let 1-L denote a sepa- 
rable Hilbert space with an inner product (. , .), and C(7~L) the set of compact 
operators acting on T~L. Then the Schatten-von Neumann classes of compact 
operators is given by 

J P (H) = {Ke C(H) : \\K\\v < oo, Vp > 1}, 

where \\K\\p = ti(\K\P) with \K\ = (K^K) 1 ^ 2 . Let A G J P {H), then for all 
z G C, the Pelmelj-Smithies formula for the p-modified Fredholm determinants 
is given by [I5j [26] 



where 



det p (l + zK) = J2^\K), 

n=0 

fp{ p) n-l • 



■ detf 



o \ 



(p) 

CP) 



y n-2 
» 



1 

» 



with 



(p) 



tr(iP) j>p 
[0 j<p-l. 

Moreover, observe that the coefficients a^\A) satisfy [27] for all p > 1 

-j n—l 

c£*\K) = - {K) V <*1 Jt cft\K) = 1. 

3=0 



(2.1) 



(2.2) 



(2.3) 



(2.4) 
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It is shown in [26 , 15] that the p-modified Fredholm determinants det p (l + zK) 
are entire functions of z whose order of zeros correspond in algebraic multiplici- 
ties to eigenvalues of the compact operator K. From Lidskii's theorem [26], [T5], 
we have for any operator K G J\{T-L) 



tr 



where A n are eigenvalues of the operator K with possibly accumulation point 
at zero, and 



tr(K) = ^2(u n ,Ku n ) 

n=l 

for any orthonormal basis {u n } n >i of 7~L. This implies that for any K G j7p(H), 
we have 

oo 
n=l 

since G Ji(H) for all j > p (cf. [15, Chap IV. 11]). There are several 
constructions for the p-modified Fredholm determinants but the one we use 
along with ( |2.1[ ) is given by (cf. [26] ) 



det p (l + zK) = Yl 



n=l 



fp-1 



(l + zA n )exp ( ^2z j (-X n ) j /j 



(2.5) 



Suppose that K G J p -i(U) then it follows from ([2j5]) that (cf. [T51I26] ) 

det p (l + zK) = det p _i(l + zK) exp ((-z) p - 1 tT(K p - 1 )/{p - 1)) . (2.6) 
Thus if K G Ji(H) then from we have (cf. [26] [25] [15] ) 

det 2 (l + zK) = deti(l + ziT)exp(-z tr(K)) . 

By a trace class operator, we refer to any operator K G J\ (H), and by a Hilbert- 
Schmidt any operator K G J2^H)- Let X C R, and consider the integral 
operator K defined on % by 



Ku(x) = / k(x,y)u(y)dx, xGl, 
Jx 



(2.7) 



where fe(a;,2/) G L 2 (X x X). Then K is compact [HI [27], since \\K\\ 2 = 
||^||l 2 (XxX) < 00 • Unfortunately, such a characterisation for trace class op- 
erator is not available. Nevertheless, if the integral operator K is of trace class 
and induced by a continuous kernel then its trace satisfies (cf. [26 ] 125] [T5]) 



tz{K) 



k{x, x)dx. 



(2.8) 



For such an operator K, its determinant deti(l + zK) is called the Fredholm 
determinant. The Fredholm determinant, as introduced originally in [14 , is an 
entire function of z which characterise the solvability of the following equation 



u(x) + zKu(x) = f{x), x G X, 
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where / and the integral kernel k associated with K are both assumed to be 
continuous in X. The latter is given bv [24l 126} H5 | [27] 



d(z) = ^z n a n (K), 



n=0 

where 



®n(K) = . I det(k(x p , x q ))p iq=1 dxi • • • dx n . (2.9) 



Note that (2.9) makes perfectly sense for any continuous kernel independently 
of whether K is of trace class or not. 



3 Convergence analysis 

Let X = [a, b] and consider the following eigenvalue problem 

(I n + zK)u{x) = 0, xe[a,b], (3.1) 



where z G C and if, defined by equation (2.7), is an integral operator acting 
on H. Basically to solve problem (3.1), there are two preferred methods; the 
expansion or projection and the quadrature method (Nystrom-type method). 
But for the sake of simplicity, we rather choose the latter method since its 
implementation is straightforward and it is efficient for smooth kernels (cf. [8]). 
However, the results in this paper are also applicable for the projection method. 
This is under the hypotheses stated in [28 which roughly assumes that K is 
compact and that both \\K — P n K || and \\K n — P n K\\ tend to zero as n — » oo, 
where P n is a projection operator. The type of kernel function k(x,y) that 
we consider in this paper are continuous everywhere in the domain except on 
the diagonal, i.e. the set {(x,y) G [a, b] x [a, b] : x = y}, and they satisfy the 
following hypotheses. 

Hypothesis 3.1. Write k x (y) = k y (x) = k(x,y). Assume that for all x G [a, b], 

k x G L^a.b], 

sup || fez || I,* < oo 
xE[a,b] 

and for x\,xi G [a, b] 

\\k Xl - k X2 \\ L i 0, as xi x 2 . (3.2) 



Observe that, we have for xi,x 2 G [a, b] and using (3.2) that 

\Ku{x{) - Ku(x 2 )\ < \\u\\ n \\k Xl - k X2 \\ L i 

Hence, it follows that the operator K maps T~L into the space of continuous 
functions, C[a,b]. Let B = {u G T~L : \\u\\u < 1}- Under Hypothesis 3.1, the 
integral operator K is compact. Indeed, the set K(B) C C[a, b] is bounded in 
T~L since for all u G 5 

\\Ku\\ n < oo, 

which follows from finiteness of the domain and the fact that 

\Ku(x)\ < sup \\k x H^i, for all x G [a, 6] and u e B. 

xE[a,b] 
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Moreover, for xi,X2 G [a, b] 

\Ku(xi) - Ku(x 2 )\ < \\k Xl - k X2 \\ L i 0, as x x -> x 2 . 

The Arzela-Ascoli lemma tells us that a totally bounded set in C[a, b] is a 
bounded equicontinuous family of functions. Hence, the set K(B) is totally 
bounded which implies that K is compact. Henceforth, we assume that the 
kernel function k(x,y) satisfying Hypothesis 3.1 is given by 

H x , y) = g( x , y)K x , y), ( 3 - 3 ) 



where g(x,y) satisfies Hypothesis |3.1| and h(x, y) is continuous everywhere in 
the domain. Note that with Hypothesis |3.l[ we can consider g{x, y) — \x — y\~ a 
with < a < 1. Assume for fixed x n G [a, 6], u G C[a, b] and N sufficiently 
large that 

N 

h(x n ,y)u(y) « [h(x n ,y)u(y)] N = ^a nj Pj(y), (3.4) 

3 = 1 

where the Pj are a polynomial basis for C[a, 6], preferably orthogonal. Then we 
define the operator Kn by 

N 

K N u(x) = ^2w j (x)h(x,y j )u(y j ), (3.5) 
where for fixed x n G [a, b] 

N 

MXn) = Y.^n) (Pjiyi))' 1 , Z = 1, . . . , JV 

i=i 

with 

PjM= / g(%n,y)Pj(y)dy. 

J a 

The operator i<0v maps H into C[a, 6], as well. This is because for all x G [a, 6], 
the continuity of Wj follows from that of (3j since by assumption the function 
g satisfies equation (3.2). Moreover, by assumption h is continuous hence the 
product Wj(x)h(x,yj) is continuous. Thus for all u G B and x G [a, 6], we have 

N 



\K N u(x)\ < ^2\wj(x)h(x,yj)\ 



iV 



< max 



ic€[a,6] . 

J = l 

< sup H^Uli max \h(x,yj)\ < oo, ViV > 1. (3.6) 

xe[a,b] xe[a,b] 



Since if h(x,y) = 1 and u(x) = 1 together with the assumption (3.4), we have 
as N goes to oo 



N 



xe[aM~[ xe[a,b] 
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Remark 3.1. A quadrature method as described in [5 can be applied in order 
to approximate h(xi,y)u(y) instead of the projection method (3.4). 



Given the operator K^, equation (3.1) is approximated by 

N 

un(x) = zKnun(x) = z^^Wj(x)h(x,yj)uN(yj), xE[a,b]. (3.7) 
i=i 

The Nystrom-type method (cf. [22 j) is then obtained by substituting x = yi in 
equation (3.7). This yields a finite dimensional eigenvalue problem given by 



N 



u N (Vi) = zK N u N (yj) = z^w j {y i )h{y i ,y j )u N {y j ), i = l,...,iV. (3.8) 

3 = 1 

To illustrate the use of the p- modified Fredholm determinants, suppose that 



K is of trace class and g(x,y) = 1 for all x,y G [a, &], in (3.3). In that case, 
Wj(x) = Wj is constant for all x G [a, 6] in (3.7). From equation (3.8), one can 
then deduce that eigenvalues of the operator Kn are precisely the zeros of the 
function (In(z) defined by 



d N (z) = det CN (I N - zK N ). 

This follows from the fact that 

deto(A - AJ n ) = (-A) n deto(/ n - zA), 

for a linear operator A G C nxn and A ^ 0, where z = A -1 , 
given in [27[ Chap VI], we can write 



(3.9) 



det C Ar(/Ar - zK N ) 



oo 

E 

n=0 



(-l) n a nN z n , 



From arguments 



(3.10) 



where a n N is defined as in (2.2) and (2.3) with tr(if n ) replaced by tr^jv^jy-). 
In fact, as mentioned in [8 and proved in [21], the series (3.10) must terminate 
at n = N. Thus, under the assumption that K is of trace class associated with 
a continuous kernel ft, we may write 



N 



d N(z) = Y[ (1 - z\ nN (K N )) 



(3.11) 



where A n Ar(i^Ar) are eigenvalues of the operator K^, for n = 1, • • • , N. Since by 
assumption the kernel ft is continuous everywhere in the domain then tr C iv (if at) 



converges to a finite tr(K) which is defined by equation (2.8). Hence for bounded 
z G C, the determinant given in (3.11 ) or (3.10) converges uniformly to the Fred- 



holm determinant d(z) (cf. [HI Theorem 6.1]). It then follows that eigenvalues of 
Kn converge to that of K. To generalise Theorem 6.1 of [8j other than for con- 
tinuous kernel, we assume that the operator for all N > 1 satisfy Anselone's 
hypotheses of collectively compact operator [lj. That is, a set {K n } n >i is called 
a collectively compact if 



Al. K and K n are linear operators on the Banach space B into itself. 
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A2. K n u — >• Ku as n — » oo, for all ^ G S and n > 1. 

A3. The set {if n }n>i is collectively compact that is, {K n u : n > 1, ||u|| < 1} 
has compact closure in S. 

Let C(H) denote the set of bounded linear operators on H. 

Definition 3.1 (Anselone pQ). A set JC C C(H) is collectively compact provided 
that the set 

KB = {Ku : K e%, u e B} 

is relatively compact^ A sequence of operators in C{T-L) is collectively compact 
whenever the corresponding set is. 



Under Hypostesis 



3.1 



the set {Kn}n>i is collectively compact. Indeed, the 



operator Kn defined in ( |3.5| ) satisfies Al. Regarding A2, observe that 



\Ku(x) — Knu(x 



)\<M max / \h(x,y)u(y) - [h(x,y)u(y)] N \dy, 

xe[a,b] J a 



where M = sup^^^ 6 ] ||^|| L i . Given the assumption (3.4), we have for suffi- 
ciently large TV 

Kjyu — >• Ku. 
For w G 5, we have by equation ( |3.6| ) that 

\\K N u\\ < M max \h(x,yj)\ < oo, V7V > 1. 

xE[a,b] 

This shows the uniform boundedness of {K^} with 

\\K N \\ < M max \h(x,yj)\ < oo, V7V > 1. 

xG[a,6] 

Furthermore, for Xi,X2 £ [a, b], we have 

AT 

|Kiv^(^i) - K N u(x 2 )\ < ^2\w j (xi)h(x 1 ,y j ) - Wj(x 2 )h(x2,yj)\ 

AT 

^ ^(iK-(^i) -^( X 2))^(^l,%)l 

|00i,%-) - /l(x2,2/j))^j(x 2 )|) 

+ IKII \H x i>yj) - h(x2,yj)\) 



iV 



Since the functions iUj and are continuous, it then follows that for all TV > 1 
\Knu(x\) — Knu(x2)\ — » as x\ — >• £2- 



1 Anselone pQ: Relatively compact, sequentially compact and totally bounded are equivalent 
in a complete space 



Hence, by the Arzela-Ascoli lemma again, the set {Kn(B)}n>i is relatively 



compact. Therefore from Definition 3.1 , the set {Kn}n>i is collectively com- 
pact. Having proved the compactness property of the operators K and {Kn}n>i, 
we are now ready to generalise Theorem 6.1 of [8 to integral operators belonging 
to Jp{T-L) for p > 1. Under the assumption that {Kn}n>i, obtained whether 
by quadrature or by projection method, is collectively compact then for N suf- 
ficiently large each eigenvalue A of K is a limit of a sequence of eigenvalues A at 
of Kn (cf. [2, 28 J. It then follows that the kth power of A is also a limit of 
the kth power of the sequence A at. Hence for bounded z, one gets a uniform 
convergence since the error in approximating the p-modified Fredholm deter- 
minants depend on the error in approximating each eigenvalue A of K which 
is uniform depending on the approximation error Ku — Kjyu. However, there 
exists some special type of integral operators K that do not fully satisfy Hy- 
pothesis [3j] but are in J±(H) (cf. [7 ). Consider for instance, operators with 
kernel k(x,y) = H(x, y)\x — y\~ a where H(x,y) is assumed to be continuous and 
a = 3/2 (cf. [7]). In that case, the convergence analysis of this paper no longer 
applies since Wj(x) is not Lebesgue integrable. However one strategy that we 
could use, as long as we are only interested in eigenvalues, is to compute the 
Fredholm determinant of the pth power of the operator K. This because for 
some the pth power integral operaor K p is of trace class and its kernel is con- 
tinuous. Consequently, it may happen that in some case a relationship between 
the Fredholm determinant corresponding to K p and the p-modified Fredholm 
determinants associated with the operator K could be established. This is of 
course will be possible depending on the the properties of det p (i^ + zK). In 
particular, if K e J 2 {U) or K G Js(H)(K e Ji{U)) we have the following 
theorem. 

Theorem 3.2. Suppose that K G J 2 (H) then 

detain - z 2 K 2 ) = det 2 (I n - zK) det 2 (I n + zK). 

Moreover, if det2(In ~ zK) = det2(In + zK) then 

(det 2 (I n - zK)) 2 = detain - z 2 K 2 ). 

If K e J 3 (H) then K 2 e J 2 (H) and we have 

det 2 (I n - z 2 K 2 ) = det 3 (I n - zK) det 3 (/« + zK) 

Proof. Since the product of two Hilbert-Schmidt operators is of trace class, we 
have 

oo 

detain -z 2 K 2 ) = H(l-z 2 \ 2 n ) 

n=l 

CO 

= Y[ ((! - z\ n ) exp(A n z)(l + z\ n ) exp(-A n z)) 

n=l 

CO CO 

= II ((1 - ZXfl ^ ex P( A n^)) Yi + zXn ^ ex P(-^n^)) 
n—1 n—1 

= det 2 (I n - zK) det 2 (I n + zK) 
= (det 2 (I n -zK)) 2 



8 



Now if K e Jz(U) then observe that K 2 satisfies (cf. [HI Theorem 11.2, Chap 
IV]) 

tr(|if 2 | 2 ) < tr(\K\ 4 ) < (tr(|K| 3 )) 4 / 3 < oo. (3.12) 
Hence K 2 is Hilbert-Schmidt. It then follows that 

oo 

det 2 (I n - z 2 K 2 ) = JJ (1 - z 2 \ 2 n ) exp(z 2 X 2 n ) 

n=l 

oo 

= I] ((1 - z\ n ) exp(A n z + z 2 \ 2 n /2) 

n=l 

x (1 + z\ n ) exp(-A n z + z 2 \ 2 n /2)) 



= detain - zK) detain + 



□ 



Remark 3.2. If if G Ja{U) then from fl3.12| ), we have 

det 2 (/H - ^ 2 i^ 2 ) = det 4 (/H - det 4 (I^ + 

And if det^i-^ — 2; K) — det^i-^ + zif ) then we have 

det 2 (/ w - z 2 K 2 ) = (deU(I H - zK)) 2 . 

More generally for if G J2p(T~L), we have for Re(z) > and p = 3, 4, . . . 

det 2p (i'H - ^f) clet 2p (i^ + K ) = c p (z) det rp/2 ] (I n - z 4 K 4 ), 

where c p (z) = 1/ det p (i^ + z 2 K 2 ). The above equality holds also for K G 
J2p-i{7~L). 

In what follows, we define for all p > 1, d p (z) := det p (l + zK) and 

00 

dNp(z) '■= det pC N(I N + ^ifi\r) = ^a^(K N )z k . 

fc=0 



Theorem 3.3. Suppose that K given in (2.7) is in J P (H) and {Kn}n>i is 
collectively compact with Kn defined by ( |3.5[ ). Then for p > 1, we have 

dN P (z) -> d p (z) 

converges uniformly for all bounded z as N — > 00. 

Proof. Let z G C be bounded by M > 0. Then we have 

00 

- d Np (z)\ < J2\ a J P) ( R ) ~ a j P N(K N )\ MK 
3=0 

Observe that from (2.4) and (2.3), the error ) — a^(ifjv)| can be deduced 

for all j > p from 



|tr(^')-tr c -(^P'- 



TV 



< 



E^-E AJ , 

n=l n=l 

AT AT 

E^-E^ 



N 



N 



E ^ 

n=AT+l 
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where A n and X n N are eigenvalues of the operators K and if at, respectively. 
As N goes to infinity, we have 



E i*»K->o- 



n=N+l 



Therefore 



E k 



n=N+l 



< E |AnK">0. 

n=N+l 



On the other hand, since {Kn}n>i is collectively compact, we have for some N 
sufficiently large, e > and for each n G {1, . . . , N} (cf. [2, 28] and [TJ Theorem 
4.8]) 

I A n — A u at| < e. 



Observe that 



l A n - Kn\ = l A n - A nA r| 



j-1 
ra=0 

< | A n - X nN \ ^2 



<j\K\ j e. 

This is because for some N sufficiently large, we have 



j'-i 
^ A^A^ 

ra=0 



TV 



: i|A„P'. 



(3.13) 
(3.14) 

(3.15) 



It then follows from (3.14) and the continuous embedding of J P (1~L) C J q (H) 
for p < q that 

E|A^-A^|<eiE|A n K< e j||if||^. 

n=l n=l 

Thus for some N large and e chosen arbitrarily small, we have 

\ti(K')-tr c ^(K 3 N )\ <e. 
Consequently, we have — d/vp(2)| < e as A/" — » oo. □ 

Remark 3.3. Suppose that fc(x, is continuous on [a, 6] 2 then (cf. [19]) we have 

tr(if n ) = / k n (x,x)dx, (3.16) 
where k\{x, y) = fc(x, y) for all x,y E [a, 6] and for n > 2 

/c n (x,:r) = / /c(x, x n _i)fe n _i (x n _i , x)dx n _i . 
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Note that in the process of computing the finite dimensional determinant (3.10), 
we essentially approximate the multiple integrals (3.16) by a product quadrature 
rule Q ni that is 



l {tr(K n ))=tr CN (K n N ) 



N 

ill- :3n- 



w 



31 



Therefore the error becomes 

\tx{K n )-tv CN {K n N ) 



\tv(K n )-Q n (tv(K n ))\, 



(3.17) 



which could be related to the quadrature error \Ku{x) — K^u{x)\ for x G [a, b]. 
Hence, if we assume that the kernel k{x,y) is smooth then one might expect an 
exponential convergence of ( |3.17| ) as mentioned in [8]. 

Under the assumption that {Kn}n>i is a collectively compact operator and 
exploiting the results of [3j [2j [28] , we can now estimate the rate of convergence in 
evaluating the determinant djsr p {z) at z = A -1 , where A G cr(K), the spectrum 
of K. In fact, this estimate depends essentially on the error in approximating 
the integral operator K in (2.7). In what follows, c is always a positive constant. 

Theorem 3.4. Assume that {Kn}n>i is collectively compact. Let A no ^ 
and X no N denote the eigenvalues of K G J V {T~L) and K^, respectively. Then for 
some N sufficiently large and for p > 1, we have 

\d p (z no ) - d Np (z no )\ < cmax{||i^ - K N Ui\\ 1/l : 1 < i < m}, (3.18) 

where z no — A" 1 , {ui, • • • ,^ m } is a basis for Ker(A no — K) 1 and m and I are 
the multiplicity and the index of \ no , respectively. 

Proof. For simplicity, we consider the case p = 1. However the following proof 
holds for p > 2, where we make use of the definition of the determinant given 
in (2.5) for finite dimensional matrices. Evaluating the Fredholm determinant 
at its zero z m = A"*, for some 1 < no < N and A no ^ 0, we have 



\d(z no ) - d N (z no )\ = \d N (z no ) 

N 

= no- 



A no \ n N 



n=l 
N 



H [Ko (^o - Kn)] 
N 



(3.19) 



Since {Kn}n>i is a collectively compact operator, it then follows from [3] and 
[28j Theorem 3] when n = no in the right-hand side of (3.19) that 



\d(z no ) - d N (z no )\ < cmax{||if^ 
for some N sufficiently large. 



K N ui 



1 1/2 • 



: 1 < i < m}, 



□ 
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For the projection method, it suffices to replace the right-hand side of (3.18) 
by the error bound in [28, Theorem 3]. 

Remark 3.4. Under normalising the eigenfunction the set {u\, • • • , u m } C B. 
On the other hand, since the set {Kn (B) }n>i and K(B) are totally bounded 
and mapping H into C[a, 6], then, K^Ui converge uniformly to Kui [TJ Propo- 
sition 1.7], that is as N — >> oo 

\\Kui-K N Ui\\ ->0. 



Theorem 3.5. Suppose K and {Kn}n>i are as in Theorem 3.4 ■ Then for 
some N sufficiently large and z, an element of the resolvent set p(K), we have 

\d p (z) - d Np (z)\ < c/3(N,l no )\\K\\P, (3.20) 

where 



f3(N,k 



m&x{\\Kui - K N Ui\\ 1/ln o :l<i< m„J 



(3.21) 



and l no = max{/ n : n = 1, . . . , N}. 



Proof Suppose that z G p(K). Then combining equation (3.13) and (3.15), and 
the error estimate of [3] for each A n , n = 1, • • • , TV, we have, for all k > p and 
for some N sufficiently large, 



N 



\tr(K k ) - tr C N(K k N )\ < ck ^|A n | fc max{||K^ - K N Uifl l « :l<i< m n }. 



Given l n , we have for all n = 1, . . . , N 



mBx{\\Kui - K NUi \\ 1/l - :l<i<m n }< f3(N, Z no ), 



with (3(N,l no ) given in (3.21). It then follows for sufficiently large N that 



\tr(K k ) - tr CN (K k N )\ < ckp(N, l no )\\K\f p . 



Hence (3.20) follows. 



□ 



Observe from equation (3.20) and equation (3.18) that for a given eigenvalue 
Zq 1 of index I <l no , we have for some N sufficiently large and z G p(K), that 



\d p (z ) - d Np (z )\ < \d p (z) - d Np (z)\. 



(3.22) 



Thus, we note from (3.18) and ( |3.20[ ) that the rate of convergence is controlled 
by the index of the eigenvalues. 

Remark 3.5. Given the above inequality (3.22), we can observe that numerical 
computation of the p-modified Fredholm determinants is similar to the numerical 
interpolation in which the interpolation points are the eigenvalues z~ x of K. 

Corollary 1. Let the integral operator K G Jp(H) and the set {ifjv}Ar>i is 
a collectively compact operator. Assume that eigenvalues X n of the operator 
K are simple, for all n > 1, and that the associated eigenvector u\ satisfies 
lm(ui) = (Re(^i) = 0). Then for some N sufficiently large and for all 
z G C and c > 0, we have 



\d p (z) - d Np (z)\ < c\\Kui - K N ui\\. 



(3.23) 
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Indeed, since the eigenvalues A n of K are simple then their algebraic multi- 
plicities m n as well as their indexes l n are equal to 1, for all n > 1 [21]. It then 
follows that 

= 11^1-^1 1|, 

and also from ( 3.22| ) that 



\d p (z) - d Np (z) \ = c\d p (z n ) - d Np (z n )\ 

where z~ x £ Hence, equation ( |3.23| ) follows. In other words, Corol- 

lary [T] tells us that the rate of convergence in computing the determinants is 
proportionally the same for any z £ C, that is for z~ x £ cr(K) or z £ p{K) 
(cf. Figure [I]). 

Remark 3.6. Suppose for example that K in Corollary [I] is associated with a 
kernel which has a jump discontinuity on the diagonal in the first derivative. 
Then for a given quadrature method which does ignore the discontinuity region, 
the rate of convergence is the same for all z £ C like in Corollary [T] However 
if the method takes into account the discontinuity then we might expect better 
convergence for z^ 1 £ cr(K) than for z £ p(K) (cf. Figure [I]). 

4 Numerical Results 

In this section we numerically evaluate the p-modified Fredholm determinants 



associated with the operator K given by (2.7), where its corresponding kernel 



k(x,y)=\ kl ^ y) : a ^ V ^ X (4.1) 
\k 2 (x,y), x<y<b 

is whether discontinuous along the diagonal or just continuous, and satisfies 
Hypothesis |3.1| For the numerical approach, we essentially follow the method 
described in [20 . To this end, let us now give a brief descprition of the method 
but more details are found in [20]. We assume that k\ and k 2 can be approxi- 
mated by Chebyshev polynomials T n (x), that is 

N 



ki(x m ,y) =^2a mn T n (y) 

n=l 
N 

feOm, Jj) = ^ CL mn T n {y), 



n=l 



for a fixed x m £ [—1,1] and m = 1, • • • , N. Then an approximation of the 
eigenvalue problem (|3.1|) is given by (cf. [20] ) 



(Vat + z b -^ (CSrC' 1 + CSiC' 1 o K 2 ) S j u = 0, (4.3) 

where o denotes the pointwise multiplication, C = T n (x m ), K\ — ki(x m ,x n ), 
K 2 = k 2 (x m ,x n ), S r and Si are right and left spectral integration matrix, 
respectively, and u = (u(xi), . . . , u(xn)) t • Hence 

dNp(z) = detAr(/Ar + zK N ), 
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where 

N 

K N u(x) = ^ (wjki(x,yj) + a j k 2 (x,y j ))u(y j ) 

with {Kn}n>1i a collectively compact operator. One of the advantages of using 
Chebyshev polynomials lie on the fact that the coefficients in the expansion of 
an indefinite integral can be easily obtain from that of the series expansion of 
the integrand, in terms of the Chebyshev polynomial [10]. Therefore with this 
property, the Chebyshev polynomials are extremely useful for integral equa- 
tions associated with discontinuous kernel function along the diagonal. In all 
the examples below, we use the Nystrom-Clenshaw-Curtis and the Nystrom- 
Gauss-Legendre referred to NCC as in [20] and NGL, respectively. The latter 
method does not take into account the discontinuity of the kernel function k 
however the NCC does. 

Example 1. For our first example, we consider the problem studied in |8], i.e 

u(x) = zKu(x) = / k(x, y)u(y)dy, xG [0,1], nGH = £ 2 (0, 1), 
Jo 

where 

x(l-y), x<y 



k{x,y) 



y(l-x), y<x. 



Our aim for this problem is not to compute the Fredholm determinant but to 
emphasize Corollary [I] and Remark |3.6| The operator for this example, 
is compact, self-adjoint and the corresponding set of operators {i^Ar} n >i is 
collectively compact since k is continuous (cf. pQ). Observe that eigenvalues of 
K, given by 

7T z n z 

are simple. Accordingly, Corollary [T] tells us that the rate of convergence is 
proportionally the same for all nonzero z, and it is of order 0(N~ 2 ) (cf. [8]). To 
confirm this, we display in Fig [l] the error \d(z) — (In(z) \ computed by the NGL 
and the NCC method at the eigenvalue, z\ = A^ 1 = 7r 2 and at z = 1. Recall 
that at the roots z n , the p- modified Fredholm determinants d p (z n ) = so the 
error is just \dN p (z n )\. 

Example 2. Our aim for the second example is to show the different rates of 



convergence shown in (3.22), for z £ p(K) and z 1 G &(K). The kernel function 



associated with the integral operator K is given by 

K X ,V) = ^-\\ X ~V\ + ^(x-y) 2 , x,y G [0,1]. 
The operator K is self-adjoint with eigenvalues given by (cf. [19]) 

= X n = ^ n 2 n 2 ' U = 1 ' 2 ' " " ' ' 

Since ^^Lil^nl < oo then the operator K is of trace class and its corresponding 
Fredholm determinant is 



^)=n(i-^) 



2 - 2cos(v^) 
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Figure 1: log 10 of the error evaluated at z\ — \ x 1 = 7r 2 (left) and at z = 1 
(right) of Example [I] 

Each eigenvalue A n has an algebraic multiplicity m n = 2, and its corresponding 
index l n = 1 for all n > 1, since if is self-adjoint. Therefore following the same 
line of arguments in [8] which led to an error of order 0(N~ 2 ) in Example [I] 
we can conclude that 

P(N, 1) = max{\\K Ul - K N ml \\Ku 2 - K N u 2 \\} 
= 0(N~ 2 ) 

where u\{x) = y/2 cos(2n7rx) and U2(x) = y/2 sin(2n7rx) (cf. [E]). However, 
numerical observation of the NGL method shows that the rate of convergence 




10° 10 1 10 2 10 3 10 4 10° 10 1 io 2 io 3 io 4 



M N 

Figure 2: log 10 of the error evaluated at z\ = A^ 1 = 47r 2 (left) and at z — 1 
(right) for Example [2] 

at the eigenvalue z\ = 47r 2 is 0(N~ A ) (see Figure (left)), whereas at z = 1 
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it is proportional to 0(N 2 ) which is also the rate for the NCC method (see 
Figure [2] (right)). 

Example 3. Let consider 

u{x) -zj vMv)dy = o, x e [-1, i\,ueH = £ 2 (-i, 1), (4.4) 



where 



k(x,y) 



1, y < x 
-1, x<y. 



(4.5) 



Note that k(x,y) is of the form given in (3.3) where g{x,y) = sgn(x — y) and 
h(x, y) = 1. The integral operator K associated with the kernel function (4.5) is 

Hilbert-Schmidt operator since 1 1 1 1 ^ 2 (( < an d it is a normal operator 

as well, that is K*K = KK* . The coefficients o4 \K) in the expression of the 
Hilbert-Schmidt determinant d,2(z) are easily computed via (2.9) by setting to 
zero the diagonal elements (cf.[T5]), and are given explicitly, for k = 0, 1, 2, ... , 

by 

1 IV, n = 2k 
n\ jo, n = 2fc + 1. 



a( 2 )(if) = 



Hence the Hilbert-Schmidt determinant is 

d 2 (z) = det 2 (/^ - zlf) = cosh(2z) 
with purely imaginary zeros satisfying 

Zn 1 = X n = ~^- 1 : , 71 = 0,1,2,... . 

7r 2n + 1 

The eigenfunctions ?i n associated with the eigenvalues A n are 



u n {x) = exp 



(i|(2n+l): 



From (4.6), it is clear that the trace of K is divergent but 

1 



^ 2 ) = -4t£ 



7T- " (2n+l) 2 



= -4. 



(4.6) 



(4.7) 



(4.8) 



The factor of 2 before the sum in ( |4.8| ) takes into account the fact that each 
eigenvalue of the operator K must be counted twice since ^(^n) = d>2(z n ) = 0. 
Observe that (4.8) can also be written as (cf. [27|[T9]) 

ti(K 2 ) = -4 = J J k(x,y)k(y,x)dxdy. 



From (2.5), we have 



d 2 (z) = n 



n=0 



(1 



4z 

V(2n + 1) 



) exp(z 



4z , 
7r(2n + 1)' 



cosh(2z). (4.9) 
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As a numerical proof of Theorem |3.3[ we plot in Figure [3] the approximate 
determinant c?2iv(^) and the original d,2(z) for \z\ < 1. In Figure [3] (left), we plot 
the imaginary part of the Hilbert-Schmidt determinant and in Figure [3] (right) 
the real part. One obviously note that as TV gets bigger, g^aK^) converge to d,2(z) 
at the rate proportional to 0(N~ 1 ) due to the discontinuity in k(x,y). Again 
to confirm that ( 3.22| ) hold, we plot the error \d,2(z) — c12n(z)\ at z = 1 G p{K) 
and zq = Aq 1 = Z7r/4. Since K is a normal operator then the index of the 
eigenvalues l n = 1 for all n > 1. Although, the eigenvalues of if are simple 
but as the eigenfunctions (4.7) are complex then Corollary [l] does not hold 
for this example. Numerically, we observe that the NGL method (dotted line) 
produces a rate of convergence proportional to 0(N~ 2 ) Figure H (left) and 
0(N~ 1 ) Figure |4] (right), for zq = iir/4 and z = 1, respectively. 




Figure 4: log 10 of the error evaluating at z\ = \ x = iir/A (left) and at z = 1 
(right) for Example [3j 
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To implement the strategy mentioned in the paragraph preceding Theorem |3.2[ 
we compute the iterated integral operator. Observe that the 2-iterated kernel 



k 2 (x, y) of (4.5) is continuous Hermitian everywhere for (x, y) G [—1, 1] x [—1, 1] 
and it is given by 

k 2 (x,y) = -2 + 2\x-y\. 

The integral operator K 2 associated with kernel k 2 is of trace class since k 2 
satisfies condition of [15j Theorem 8.2, Chap IV] with a > 1/2. Therefore its 
trace is 



tr(K 2 ) = J k 2 (x,x)dx = -4, 



and its Fredholm determinant is 



d(z) = det(I n - z 2 K 2 ) = i(cosh(4z) + 1). (4.10) 



Note that d 2 (z) satisfy the conditions given in Theorem 3.2 Hence we must 
have that 

(d 2 (z)) 2 = d(z). (4.11) 
Indeed the above equality holds since 

(cosh(2z)) 2 = i(cosh(4z) + 1). 



Note that due to the discontinuity in the kernel function (4.5), the approximate 



Hilbert-Schmidt determinant d 2 ]y(z) converges slowly to d 2 (z). However given 



the relation (4.11) and that k 2 is continuous, one can improve the order of 
convergence in approximating d 2 (z) by approximating det(i^ — z 2 K 2 ) instead. 
Indeed, k 2 has similar properties as the kernel in Example [2] Therefore, the 
same behaviour of the error can be expected as shown in Figure [2] 

Remark 4.1. The same behaviour of the error holds for the self- adjoint operator 
iK with eigenvalues ±iA n . 



Example 4. We consider equation (4.4) with kernel given by 

k{x, y) 



\x-y\<*' 



The kernel function k(x,y) is of the form given in (3.3) with h(x,y) = 1. Hence 
the integral operator K associated with the above kernel is compact, self-adjoint 
[2] and positive. In particular, for < a < 1/2 the integral operator K is 
Hilbert-Schmidt (cf.p~9[ |6]). Therefore, the 2-modified Fredholm determinant 
is computed as in Example [3] by substituting zero in the kernel function for 
x = y [E]. Then we shall focus on the case a = 1/2. For this case, the 
integral operator K is not Hilbert-Schmidt since ||&||l 2 ((-i,i) 2 ) is unbounded 
hence K G J P {T-L) for p > 3. Since K is self-adjoint and positive operator so is 
the integral operator K 3 . Moreover, the 3-iterated kernel k^(x^y) is continuous 
(cf.[6]) it then follows that the integral operator K 3 is of trace class (cf. [26], 
Theorem 2.12]) with 



tr(K 3 ) = J k 3 (x,x)dx = ||K 3 ||i < oo. 
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This implies since K is self-adjoint and positive operator that 



tr(\K\ 



(tr^ 3 )) 1 / 3 



< 00. 



Hence K is in J^{T~L). For the numerical computation of the 3-modified Fred- 
holm determinant ds(z) = dets(Iy^ — zK), we need to compute numerically 
the eigenvalues of K and form the finite dimensional version of equation (2.5). 



However, given Theorem 3.2 we are not required to, we only need to have an 
explicit expression of the 2-iterated kernel k 2 {x,y) and set its diagonal values 
to zero. Using Maple, the 2-iterated kernel k 2 (x,y) is 



k 2 (x,y) 



2v/(l-2/)(l-x)) +7T 
■ + 2^/(1 + 2/)(l + a;)), 



• In (2 - y - : 

+ In (2 + y 
■ In (2 + y + x - 2 v / (l + y)(l + x)) + - 

+ ln (2-y 



■ + 2y/(l-y)(l-xj), 



x < y 



x>y 



To confirm Theorem |3.2| for K G Js(Jl), we display in Figure [5] (left) the plot 
of dsN(z) = det3Ar(ijv — zK^) computed using the five eigenvalues of largest 
modulus and d 2 n (z) = det2Ar(/jv — z 2 K^) computed by the rectangular rule. 
In Figure [5] (right), we plot dsN(z)dsN(— z) and d 2 N(z). Since K is positive 
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Figure 5: Numerical approximation of ds(z) (solid line) and d 2 (z) (left), and 
ds(z)ds(—z) and d 2 (z) (right). 

operator then its eigenvalues are positive. Therefore from Theorem |3.2[ we 
have for real z > 

d3N(z) = c N (z)d 2N (z), 

where cn(z) = l/d3jv( — z) is a nonzero entire function. From Figure [5j we 
observe that the smallest root z\ = A^ 1 of ds n (z) is simple, and together with 
the self-adjointness property of K we conclude that eigenvalues of K are simple. 
Moreover, we have observed numerically that the eigenfunctions of K are real. 
Hence the error in evaluating the 3-modified Fredholm determinant satisfies 
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Corollary [T] that is for a fixed real z > we have 



\d 3 (z) - d 3N (z)\ = \c(z)d 2 (z) - c N d 2N (z)\ (4.12) 
KcWKm-KNmW 

where c(z) = l/d 3 (—z). Since the error is proportionally the same then we can 
bound the above error with the error in approximating the 2-iterated operator 
K 2 , that is 

\d 3 (z) - d 3N (z)\ < \c(z)\\d 2 (z) - d 2N (z)\ 

< c\c(z)\\\K 2 Ul - K^ml 

It then follows from the uniform convergence of the Hilbert-Schmidt determi- 
nant det2 that equation ( 4.12| ) convergences uniformly too. For numerical proof 



(cf. Figure[6|, we plot the error \d 2 (z) — d 2 N(z)\ using the rectangular rule eval- 
uated at z = 0.1 which according to the graph is of order 0(N~ 1 log(TV)). 




Figure 6: log 10 °f the error \d 2 (z) — d 2 N(z)\ where z — 0.1. 



5 Conclusion 

In this paper, we have given theoretical and numerical results concerning the 
approximation error (d p — djsr p ), where d p is the p- modified Fredholm deter- 
minants and d^p is the finite dimensional determinants corresponding to d p . 
These results are as follows; first, we have shown that the approximation error 
(d p — dN P ) is uniform in a bounded domain. Second, we have given the rate of 
convergence when evaluating (d p — d^p) at an eigenvalue or at a point of the 
resolvent set. As a consequence, we have observed that numerical evaluation 
of the p-modified Fredholm determinants is nothing else than an interpolation 
in which the interpolation points are the eigenvalues of the operator K. From 
the well known result of interpolation theory which states that if the interpo- 
lated function is continuous then the error is uniform in a bounded domain, this 
confirms the uniform convergence obtained in Theorem |3.3| Although we dealt 
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with a finite domain of R, an extension of the present analysis to X C W 1 is 
possible. This is of course under the hypothesis that the kernel k(x,y) is such 
that the set {Kn}n>i is collectively compact. However, numerically we believe 
that for L 2 -kernel function k{x,y) on the real line R and defined as in (4.1 ), the 
NCC method will not be effective. This is due to the unbounded nature of the 
kernel functions k\ and &2- 
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